Calculated Momentum Dependence of Zhang-Rice States in Transition Metal Oxides 
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Using a combination of local density functional theory and cluster exact diagonalization based 
dynamical mean field theory, we calculate many body electronic structures of several Mott insulating 
oxides including undoped high Tc materials. The dispersions of the lowest occupied electronic states 
are associated with the Zhang-Rice singlets in cuprates and with doublets, triplets, quadruplets and 
quintets in more general cases. Our results agree with angle resolved photoemission experiments 
including the decrease of the spectral weight of the Zhang-Rice band as it approaches k=0. 



Quasiparticle excitations in insulating transition metal 
oxides (TMOs) such as classical Mott-Hubbard systems 
or undoped high temperature superconductors (HTSCs) 
have been puzzling electronic structure theorists for 
many years [l|, ■ While photoemission experiments in 
these materials show ^] the existence of the d-states lo- 
cated both right below the Fermi energy and at much 
higher binding energies (typically ^10 eV), it is difficult 
to understand this genuine many-body redistribution of 
the spectral weight using calculations based on a static 
mean field theory 0, such as the density functional 
theory (DFT) in its local density approximation (LDA) 
0. Modern approaches, such as LDA-I-U 0, can differ- 
entiate between charge-transfer and Mott-Hubbard na- 
tures of these systems but still have difficulties in 
recovering insulating behavior of the paramagnetic (PM) 
state and tackling more complicated many-body features 
such as Zhang-Rice (ZR) singlet of HTSCs d 0. Only 
most recent developments based on a combination of lo- 
cal density approximation (LDA) and dynamical mean 
field theory (DMFT) pUi] have started to address those 



issues 



In the present work, using a novel implementation of 
LDA plus cluster exact diagonalization based DMFT we 
demonstrate how to obtain accurate spectra of transition 
metal oxides and, in particular, describe full momentum 
dependent low-energy excitations associated in those sys- 
tems with antiferromagnetic (AFM) Kondo-like coupling 
between a spin of oxygen hole injected by photoemis- 
sion process and a local magnetic moment of the tran- 
sition metal ion. These narrow energy bands are com- 
posed from the well known Zhang-Rice singlet states in 
cuprates ^ , which have recently renewed their attention 
in connection with the disappearance of their spectral 
weight as the wave vector approaches the Brillouin Zone 
(BZ) center, and the observed high energy kink entitled 
as "waterfall" feature n4| . Zhang-Rice doublets have 
been discussed in NiO [l5|, and their further general- 
izations to triplets (CoO), quadruplets (FeO) and quin- 
tets (MnO) all naturally emerge from our LDA-fDMFT 
calculations. We find that the ZR states exhibit a sim- 
ilar behavior in all systems including the loss of their 



spectral weight at the F point, which can be understood 
as the lack of hybridization between transition metal d 
states and neighboring oxygen p states, the effect most 
pronounced in HTSCs. There is a generally good agree- 
ment between our results and angle resolved photoemis- 
sion (ARPES) experiments. 

In our self-consistent LDA+DMFT calculations, the 
LDA one-electron Hamiltonians are supplemented with 
the self-energies of the cZ-electrons that are extracted 
from the cluster exact diagonalizations involving a d-shell 
of a transition metal ion hybridizing with nearest oxygen 
orbitals. We also include the effect of inter-site self- 
energies in the study of cuprates including Sr2Cu02Cl2 
and La2Cu04 by considering clusters of two copper id 
orbitals. The hybridization functions in our calculations 
are fit using a single pole approximation A{uj) = V"^ j (uj — 
Ep) which delivers the effective position of the bath level 
Ep (interpreted as the oxygen p band), as well as its hy- 
bridization V with the d states. The cluster self-energies 
are fit into a rational form E((jj) = I](oo)+Xlj Wi/(a;— Pi) 
using interpolation by three poles JJj]. Such fits allow 
to simplify the process of extracting parameters for the 
cluster Hamiltonians and performing full self-consistency 
using the LDA+DMFT for the actual materials consid- 
ered in our work. Conventional values for the Coulomb 
interaction matrices, which were found from constrained 
density functional calculations are used to perform 
the calculations at temperatures above the long-range 
magnetic order. 

Since dimensions of the many-body Hamiltonians 
quickly become prohibitively large to handle by standard 
diagonalization algorithms, we have newly implemented 
a kernel polynomial method (KPM) for extracting spec- 
tral functions in our cluster calculations Similar to 
the Lanczos method, the KPM is an iterative procedure 
which allows to recover moments of many body densities 
of energy states. Using 200 moments or so, the one- 
electron Green's functions are found to converge accu- 
rately for our cluster Hamiltonians with the dimensions 
up to 100,000. 

We first discuss a long-standing problem of low en- 
ergy excitations in HTSCs. Early LDA calculations 
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FIG. 1: (color online) Calculated low energy excitations in 
Sr2Cu02Cl2 (top) and La2Cu04 (bottom) using LDA plus 
single-site DMFT (red line) and two-site cluster DMFT (blue 
line). The experimental data are from: Ref. [l^ (open cir- 
cles), Ref. [2(^ (open triangles), Ref. i2l| (open squares). The 
blue linewidth shows the oxygen content, while the numbers 
indicate the amount of electrons in the ZR band. 



were unable to reproduce even the magnetic behavior 
of those systems [2]. Although this problem was later 
solved by the LDA-fU method [zl, these calculations 
miss the important ZR physics and are unable to ad- 
dress the low-energy scale. On the other hand, t — J and 
t~t' — t" ~ J model calculations have successfully repro- 
duced the ZR dispersion as compared to the ARPES data 

isl Til 13) 21 , 23] ■ A hot discussion has appeared in the 
literature around most recent ARPES experiments which 
have shown that the ZR band vanishes as the wave vector 
approaches the P point, with the spectral weight seem- 
ingly transferring to higher-energies as the "waterfall" 

14| . While it could be simple photoemission matrix ele- 
ment effect, most recent dynamical cluster approximation 
(DCA) based simulations of the 2D Hubbard model show 
[23| large imaginary part of the self-energy, resulting in 
asymmetric spectral functions resembling the waterfall. 
As we generally see, the HTSCs still represent a chal- 
lenge for the theory, particularly for modern electronic 
structure calculations, which try to incorporate all hop- 
ping integrals accurately and use realistic values of the 
Coulomb interaction parameters. 

In the present work LDA-I-DMFT calculations are 
performed for Sr2Cu02Cl2, whose ARPES experiments 
clearly identified the ZR band and for La2Cu04. 
First, the calculations were done for a single impurity 
case when the cluster is made of a copper ion surrounded 
by 4 oxygen orbitals. This result is presented in Fig[TJtop 
panel) by red line together with various photoemission 



experiments denoted by symbols. The dispersions of the 
ZR band in Sr2Cu02Cl2 correctly follows the ARPES 
data, although its bandwidth is overestimated by a factor 
of two, especially along (0, 0)-(7r, tt) direction. The situa- 
tion is somewhat unclear along (0, 0)-(7r, 0) line since ex- 
periments by two groups show quite different bandwidths 
(triangles [20] vs circlcs[19] and squares [2lj ) . From the 
standpoint of the single-impurity model, it is however 
clear that the dispersion of the ZR band above the Neel 
temperature is primarily governed by the short-range 
antiferromagnetic correlations in that once the hole is 
moving in the lattice it is dressed up if the short-range 
order is present. This effect is missing when the single- 
impurity approximation is adopted as the DMFT treats 
the PM regime as completely disordered state, resulting 
in the ZR band being too wide in our calculation. To in- 
clude short range magnetic correlations, we subsequently 
performed calculations with two impurities by consider- 
ing a cluster of two copper atoms surrounded by 7 oxy- 
gens (one shared orbital). The cluster exact diagonaliza- 
tion now delivers both the on-site I]ii(w) and inter-site 
I]i2(w) self-energies, which can be Fourier transformed 
to the form S(k,Ci;) = Eii(a;) -|- Si2(w)(cos fc^; -I- cosfcj,) 
and acquire the much needed k-dependence. By fixing 
the frequency w to the position of the ZR band we are 
now able to add I](k,w) to the LDA Hamiltonian. The 
result of such a calculation is depicted by the blue lines in 
FigU A remarkable band narrowing now occurs where 
the ZR bandwidths come about 0.3 eV for Sr2Cu02Cl2 
and 0.4 eV for La2Cu04. Using the calculated ZR dis- 
persion, the parameters for t — J model are also evalu- 
ated. For Sr2Cu02Cl2 they appear to be J = 0.135eV, 
t = 0.34eV, t' = -0.12eV and t" = O.OSeV, and for 
La2Cu04 the numbers are J = 0.175eV, t = 0.44eV, 
t' = -0.15eV, t" = O.lOeV. The reduction of the ZR 
bandwidths by including k-dependence of the self-energy 
demonstrates the importance of going beyond the single- 
site approximation in describing the low-energy excita- 
tions in cuprates. 

We further discuss the effect of disappearance of the 
ZR spectral weight around the P point. In FiglU the 
blue fat line on the bottom panel shows the amount of 
oxygen p states presented in the ZR band of La2Cu04. 
It is clearly seen that the spectral weight gets smoothly 
reduced as the k vector approaches the BZ center. 
The same effect is found in our calculations for other 
HTSCs. The reduction can be simply understood using 
the 3-band tight-binding argument: the k-dependent 
hybridization between dx2-y2 and Px,Py orbitals is ex- 
actly equal to zero at P, forbidding the formation of the 
ZR singlet. We can also track this effect with E(k, u>) us- 
ing a recently proposed matrix expansion algorithm [l6| : 
Assuming that two poles of the self-energy are needed to 
reproduce both the Hubbard and the ZR bands, the re- 
sulting Hamiltonian is 5 x 5 in the extended pole space (3 
for d^2_y2 ^px,py and 2 for the self-energy poles). It can 
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tively. The 3d—2p hybridization is responsible for the Bd- 
character peak just below the Fermi level [2^. The dis- 
cussion of a generalized spin-fermion model derived for a 
slab of NiO shows a well-developed set of ZR states which 
agree reasonably well within the low-energy features seen 
in the ARPES experiments 
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FIG. 2: (color online) Comparison between calculated quasi- 
particle dispersions using LDA+DMFT and the photoemis- 
sion data [a, [2^ for paramagnetic state of NiO. The blue 
linewidth and the numbers show the oxygen content and the 
amount of electrons in the ZR band. 



be simply diagonalized and the same effect is monitored. 
The amount of electrons in each band is less than one 
since auxiliary pole states are added to the Hamiltonian 
[1^ . FiglT] shows the actual number of electrons in the 
ZR band as a function of k, decreasing towards F. Here, 
all spectral weight gets transferred to the Hubbard bands 
located at the scale of U and also to the oxygen band. 

While we are unable to account for the effect of the 
imaginary self-energy in the exact diagonalization for 
small clusters, recent DCA simulations of the Hubbard 
model [23*1 have shown very asymmetric k-dependent 
spectral functions resembling the waterfall. On top of 
that a lot of other states are present below —1 eV, as 
shown in Fig[l] (bottom panel). In our opinion this pro- 
vides most natural interpretation of the recent ARPES 
data. As small hole doping occurs, the Fermi level sits 
at top of the ZR band which would again contradicts 
with the standard LDA like description of the cuprates, 
although we are unable to study the regime of optimal 
doping by the present procedure requiring too large clus- 
ters to diagonalize. 

Let us now turn to the discussion of NiO which has 
been widely chosen as another classical example of a 
strongly correlated system. Two decades ago Fujimori et 
al 24 1 and Sawatzky et al [25[ have shown that the va- 
lence band photoemission spectrum of NiO can be under- 
stood from a configurational interaction (CI) approach, 
indicating that the electronic properties here are local. 
The 3d electrons repulsion gives rise to the lower and 
upper Hubbard bands located at —9 and 4 eV, respec- 



However, NiO exhibits 
both correlation effects and band-structure effects, plus 
antiferromagnetic order under the Neel temperature of 
520 K. Below Tjv extensive studies of NiO have been per- 
formed using the LDA-I-U method, and the positions of 
the Hubbard bands have been predicted correctly when 
using the values of on-site Hubbard interaction U = 8 
eV [7|. Unfortunately, much of the d-electron spectral 
weight just below the Fermi level is found to be lacking 
in this type of calculation unless the value of U is re- 
duced to 4 ~ 5 eV making this system of Mott-Hubbard 
type rather than of charge-transfer type. Better than 
LDA+U, by using atomic (Hubbard I) approximation for 
the self-energy, one can get insulating solution even for 
paramagnetic case llj. Also, recent comparisons between 
LDA calculations of NiO in its paramagnetic (predicted 
to be metallic) state and the ARPES experiments have 
indeed revealed [5] a rather good agreement between each 
other in the directions of BZ where the LDA band dis- 
persions remain insulating. One thus faces a dilemma 
of matching different portions of the photoemission data 
with three different types of electronic structure calcula- 
tions (LDA+U with two values of U as well as straight 
LDA). 

We now show that all these problems can be overcome 
by performing LDA-I-DMFT calculations. Figl^] shows 
the comparison between calculated quasiparticle disper- 
sions and photoemission data for NiO [2^ . The Ni 3d 
spectral weight is redistributed between the lower Hub- 
bard band located at — 8 '-^ —10 eV, upper Hubbard band 
located at around 4 eV, and strong peaks located just be- 
low the Fermi level. Having the ground state in NiO as 
\p^d^), the lowest excitation here is the result of the inter- 
action between the spin of the oxygen hole and the two 
hole S' = 1 state in the d-shell. The total spin of such 
cluster can be 1/2 (doublet) or 3/2 (quadruplet). The 
antiparallel Kondo-like coupling implies that the doublet 
state is lower in energy, in accord with early predictions 
using the spin-fermion model [l5| . The situation is akin 
to the singlet and triplet situation in cuprates. The di- 
rect energy gap is seen in Fig[^ to be around 3.5 eV in 
our calculation (the experimental value is 4 eV), although 
most of the optical transitions would start at the energies 
4 eV or so as is evident from our density of states plots. 
We also see that both the positions and the dispersion 
of the oxygen bands located at the energies —4 ~ — 7 eV 
below the Fermi level agree well with the ARPES data. 

It is remarkable that similar to the HTSCs the reduc- 
tion of the spectral weight of the ZR doublet is also seen 
in NiO when the k vector approaches the F point. Here 
the same tight-binding argument applies: in the nearest 
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FIG. 3: (color online) The same as Fig. [2] but for CoO. Pho- 
toemission data are from Ref. [2?! . [28| ] 



neighbor approximation the k-dependent hybridization 
between oxygen p and orbitals vanishes prohibiting 
the creation of the ZR state. The amount of actual p 
electrons in the ZR state is shown by fattening this band 
in Fig. [21 While it decreases when k — > 0, there is some 
residual amount around F due to other hoppings which 
results in the actual number of electrons (shown by num- 
bers) to be not as drastically reduced at F as in the case 
of HTSCs. 

Our result for CoO demonstrates very similar behavior 
to that of NiO. Fig [3] compares the calculated quasiparti- 
cle dispersions with the photoemission data [13, [1^ . The 
calculated energy gap is ^2.5 eV, and the experimen- 
tal value is about 2.4 eV. Remarkably similar to NiO, 
the d-electron spectral weight is distributed between the 
Hubbard bands found at energies —9 ~ —10 eV and at 
4 eV, as well as the ZR band just below the Fermi level. 
The lowest excitation here can be described as a spin 
triplet state because of the antiferromagnetic interaction 
between the oxygen hole and the 3 hole S' = 3/2 state of 
the Co ion. We also see that the spectral weight in the 
ZR triplet band gets reduced as one moves towards the 
BZ center exactly as in the previous cases. 

We finally mention our calculations for FeO and MnO 
where a very similar behavior was found. Since the spin 
momenta on Fe and Mn sites are roughly equal to 2 and 
5/2, respectively, the corresponding generalizations of 
the ZR states are spin 3/2 (quadruplet) for FeO, and spin 
2 (quintet) for MnO. We should however mention that 
moving towards middle of the 3(i-metal oxide series, the 
hybridization effects become more and more pronounced 
and our calculations based on small cluster diagonaliza- 
tion may be less accurate. 



In conclusion, we have performed cluster exact diago- 
nalization based LDA-I-DMFT calculations for quasipar- 
ticle spectra in selected transition metal oxides. Low 
energy excitations were compared in details with the 
ARPES experiments including the ZR states in selected 
HTSCs, as well as doublet, triplet, quadruplet and quin- 
tet states in NiO, CoO, FeO, MnO, respectively. Gen- 
erally good agreement between the theory and experi- 
ments was found using single site DMFT approximation 
for classical Mott-Hubbard systems, while for the case 
of the cuprates, going beyond the single impurity is es- 
sential. The reduction of the spectral weight in the ZR 
band when approaching the F point was monitored for 
all materials and explained based on the tight -binding 
argument. 
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